Entanglement Induced Phase Transitions 
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Starting from the canonical ensemble over the space of pure quantum states, we obtain an integral 
representation for the partition function. This is used to calculate the magnetisation of a system of 
A'^ spin-i particles. The results suggest the existence of a new type of first order phase transition 
that occurs at zero temperature in the absence of spin-spin interactions. The transition arises as 
a consequence of quantum entanglement. The effects of internal interactions are analysed and the 
behaviour of the magnetic susceptibility for a small number of interacting spins is determined. 
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In classical statistical mechanics, a phase transition is 
a phenomenon characteristic of systems having internal 
interactions Indeed, if the form of the internal energy 
is reasonably idealised, there is a wide range of models 
describing phase transitions at finite temperatures that 
can be solved exactly The situation is similar for 
quantum systems, though exact solutions are typically 
difficult to obtain in this case 

One distinguishing feature of quantum systems in this 
context is the existence of entanglement. This gives rise 
to a colossal increase in the state space volume for com- 
bined systems. Furthermore quantum entanglement can 
be viewed in itself as a form of internal interaction be- 
tween the constituent particles of the system. It is nat- 
ural therefore to inquire whether entanglement has any 
role to play in critical phenomena. 

In this letter we consider an entangled system of dis- 
tinguishable spin-i particles, weakly interacting with a 
heat bath. A general integral formula for the partition 
function is obtained and used to analyse the magnetisa- 
tion and the magnetic susceptibility of the system. In 
particular we show that the magnetisation of N nonin- 
teracting particles in the low temperature limit behaves 
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where is the magnetic moment of an individual parti- 
cle, B is the applied magnetic field, T is the temperature 
of the heat bath, and ks is Boltzmann's constant. The 
A'^-dependence of the magnetisation points, in the ther- 
modynamic limit TV— >oo, to the possible existence of a 
first order phase transition at zero temperature in the 
absence of spin-spin interactions. This is a phenomenon 
with no classical counterpart, indicating that entangle- 
ment may have an effect similar to that of a mean-field 
which enhances critical phenomena. 

Let ■0" be a state vector in a complex Hilbert space TL. 
A physical observable, such as the Hamiltonian Hp, is 
represented by a linear operator acting onTi. If a system 
is in the state ip"', then the quantum expectation of an 



operator in that state is {F) = ^aF^-ip^ /-ipjip^ , where 

-ipa is the complex conjugate ofip". The expectation (F) 
is unaltered under transformations of the form tJj" — > 
Aijj" , where A is any nonvanishing complex number. In 
other words, the Hilbert space formulation of quantum 
mechanics carries an extra complex degree of freedom 
given by the overall scale and phase. It is convenient, 
therefore, to introduce the space of equivalence classes 
of vectors in Ti. modulo complex scale transformations. 
This is the projective Hilbert space VTC: a point in PH 
corresponds to all the points on a ray through the origin 
of H, except the origin itself. We refer to VTC as the 
space of pure quantum states. 

For definiteness, let us take H to be an (n + 1)- 
dimensional complex Hilbert space. Then the corre- 
sponding state space is the rt-dimensional complex pro- 
jective space CP". The state space CP", when regarded 
as a real manifold P of dimension 2n, is equipped with 
a natural symplectic structure, as well as a Riemannian 
structure known as the Fubini-Study metric Each 
point x in P represents a ray in TC. Conversely, suppose 
we assign to each point x in P a nonvanishing element 
ip^i^x) in the ray corresponding to x. Then given any ob- 
servable F acting on 7i, we can construct a real-valued 
biquadratic function F(x) on P given by the expectation 
of P^ in the state ■ip°'(x). It is straightforward to ver- 
ify that F{x) is independent of the particular choice of 
cross-section ip"{x). 

When we shift emphasis from H to P, the role of the 
complex number is taken over by the symplectic struc- 
ture, given by a nondegenerate skew-symmetric tensor 
field, and the Schrddinger equation can be expressed in 
Hamiltonian form. The formulation of quantum mechan- 
ics on the state space P thus bears a striking resemblance 
to classical Hamiltonian mechanics, with the additional 
constraint that the Hamiltonian H{x) in quantum me- 
chanics is of the special form H{x) = Hgll^(x), where 
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is the projection operator corresponding to the pure state 
xeT. 

The fact that the space F of pure states is the quan- 
tum phase space makes F the appropriate space in which 
to consider the effects of quantum entanglement. In this 
letter our interest is in the weak coupling limit between a 
quantum system and an environment with which it is in 
thermal equilibrium. Specifically, we would like to know 
how to characterise thermal equilibrium states on F, that 
is, to find an ensemble p{x) over F that describes the 
distribution of wave functions. This question has been 
addressed in |^], where it is shown that, if we assume 
the resulting distribution over F maximises the Shannon 
entropy associated with p{x), then the equilibrium en- 
semble is given by the Gibbs measure 



where f3 = l/fc^T, and 
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is the partition function. Here dV is the volume element 
arising from the Fubini-Study metric on F. More specif- 
ically, we have 



for the phase space volume element, where 
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and eaf3-y---s is the totally skew-symmetric tensor. 

The interpretation of the distribution (j^) is as follows. 
For each pure state x G F we compute the expectation 
H{x) of the Hamiltonian operator, conditioned on that 
pure state. For thermal equilibrium the probability that 
the quantum system is in the pure state x is determined 
by the Boltzmann weight (|^). Then the unconditional 
expectation of the Hamiltonian in the thermal state is 
given by the internal energy 



U{p) = j^H{x)p{x)dV, 
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which is equivalent to the trace formula H^p^, where p^ 
is the density matrix associated with the distribution (|^), 
given by 



p{x)ni{x)dv. 



It follows by a standard identity that U = —d\nZ/df3. 
Therefore, we would like to obtain an explicit expression 
for the partition function Z{j3), in order to determine 
properties of the thermodynamic functions. 



For the computation of the partition function, it is 
convenient to revert to the homogeneous coordinates ij)" 
{a = 0, 1, • • • , n) on the complex projective space CP". 
Then, we observe that the integration in (^) can be lifted 
to C""*"^ with a spherical constraint V'aV'" = !> which 
gives us 



Moreover, if we choose a basis such that the Hamiltonian 
operator is diagonal and substitute the standard identity 
5{x) = (27r)-i / e'«^d^ for the (5-function, then the C"+i- 
integration becomes a Gaussian and we obtain 
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where Ek are the energy eigenvalues. In obtaining ( [lO| ) 
we transform from complex coordinates ip"' to real co- 
ordinates, which makes the integration involved for each 
energy eigenvalue Ek in a double Gaussian. If we 
analytically continue f to the lower half-plane, then the 
integrand decreases exponentially and we can close the 
contour and apply the residue theorem to evaluate the 
integral. The result is 
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where we have written A = — i^ and we have discarded 
the physically unimportant factors of 27r. 

Formula ( pT| ) for the partition function is valid for any 
quantum system that can be modelled by a finite dimen- 
sional state space. In particular, if there is no degeneracy 
in the energy eigenvalues, then we find that the partition 
function reduces to the following simple expression: 
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This formula is applicable, for example, to the case of a 
single spin-^ particle in a magnetic field, for which we 
can write Ek — ^(-f — k)p,B, where fc = 0, • ■ • ,n. Then 
we have 
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for the partition function. Interestingly, this is identical 
to the expression one obtains for n classically indistin- 
(8) guishable, independent (disentangled) spin-i particles in 
a magnetic field. In the spin-i case (i.e. n= 1), as was 
addressed in Q , the magnetisation energy is given by 



U = ksT - -pB coth iPpB. 



(14) 
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Let us now apply the result (|Tl|) to a system of dis- 
tinguishable spin-i particles (e.g., electrons on a lattice) 
in a magnetic field. Classically, if the particles are not 
interacting, then the resulting partition function factors, 
and wc obtain the same magnetisation per particle as 
in the single particle case. However, quantum mechani- 
cally, this is no longer the case, because of the existence 
of quantum entanglement. That is, even in the absence of 
direct interactions, the presence of entanglement implies 
that the particles are not independent. For example, if 
the system consists of two spin-i particles, any entangled 
state (such as the singlet state) gives rise to a nonvanish- 
ing Boltzmann weight through formula (^). 

For N noninteracting spin-^ particles, the Hamiltonian 
operator is 

JV 

H ^-liBY^s.^^^fxBS,, (15) 

i=l 

where S = X^i^i^*- Note that, although the total 
number of eigenstates is given by 2^, the correspond- 
ing eigenvalues are highly degenerate. In particular, 
there are only A'^-f 1 distinct energy eigenvalues, given 
by Efc = — (-y — /c) /i_B, where the index k runs from to 
N . The degree of degeneracy associated with the eigen- 
value Efc is TvCfc. Using this energy spectrum, we have 
determined the magnetisation of the system per particle 
for N up to five, with the resulting plot given in Fig. 1. 

The results in Fig. 1 show that, even in the absence 
of spin-spin interactions, the existence of quantum en- 
tanglement gives rise to nontrivial thermal expectation 
values — a result that has no classical counterpart. Fur- 
thermore, it also indicates that, in the thermodynamic 
limit N ^ oo, the gradient of the magnetisation at zero 
temperature diverges. If this were the case, then we 
would obtain a first order phase transition at T = such 
that the value of magnetisation would jump from to 1/2 
as T I 0. However, because the formula (|l^) is not guar- 
anteed to be valid in the thermodynamic limit, we cannot 
prove the existence of a phase transition. Nevertheless, 
we may consider the asymptotic behaviour of the mag- 
netisation in the low temperature limit, and study how 
the resulting expression depends on TV. In particular, 
as r — > the contribution from the residue arising from 
the ground state energy dominates the partition function. 
Thus we find 

^ 1 

= apfiPfiBy-^" exp (16) 

where oat is independent of (3. Then using M — 
/3~^9 In Z/ dB, we obtain the behaviour of the magnetisa- 
tion as given in (|l|). This result is valid for an arbitrarily 
large but finite N, in the limit T — > 0, and agrees with 
the numerical results in Fig. 1. In particular, the coeffi- 
cient of T diverges exponentially in N, strongly pointing 
towards the existence of a transition. 
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FIG. 1. The magnetisation per particle for noninteracting 
spin-i particles. Classically, the result is independent of the 
number A'^ of particles. Quantum mechanically, as A'^ is in- 
creased, the effect of entanglement accumulates to give quan- 
titatively different behaviours for different A'^. 

The phase transition implied by the foregoing analysis 
is, of course, rather artificial, in the sense that the transi- 
tion takes place at zero temperature, and that the mag- 
netisation can take only two distinct values. This is not 
surprising however because the model we have considered 
has no spin-spin interactions. If we include interactions 
in order to have a finite temperature transition, then we 
expect the decrease of M/N with increasing A^ at fixed 
T, as observed in Fig. 1, is suppressed. 

We now turn to study a highly interacting system 
in which derivations of thermodynamic functions are 
analytically tractable. In particular, we consider the 
quantum Ising model on a complete graph, rather than 
the conventional square-lattice models. Therefore, the 
Hamiltonian operator is now given by 

N 

H = -iiB ^s,^ - Si ■ Sj 

1=1 i>j 

= -fiBS, - i J(^2 - f A^), (17) 

Without loss of generality, we can assume that the total 
number A^ = 2n of the spins is even. Then ( pi] ) gives us 

residues \s=Om= — s ^ ^ " / 

(18) 

where 

e„(s) = -m/iS - i J {s{s + 1) - f n) (19) 

and d2nis) = 2nCn-s ~ 2riC'„_s_i. Thc grouud state 
energy, in particular, is given by Eq — e„(7i). 
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With the expression (|18|) at hand, we can consider the 
behaviour of thermodynamic functions. In particular, 
for a system of interacting spins, we can determine the 
second moment, namely, the magnetic susceptibility of 
the system. Using the definition x — {^/N)dM/dB for 
the magnetic susceptibility, we analysed numerically the 
behaviour of x for a range of values of N . The results are 
shown in Fig. 2, indicating the vanishing of the second 
moment at T = and T — > oo, as well as a peak at finite 
temperature, the sharpness of which increases with TV. 
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FIG. 2. Magnetic susceptibility of the system of interact- 
ing spin-i particles on a complete graph for A'' — 1, ■ ■ • ,5. 
The numerical value J/^iB = 0.2 has been chosen for the 
interaction strength. 

There are two regimes of (|l^) which can be explored 
analytically. The first is the case of weak coupling in the 
low temperature limit. More precisely, when J/fiB < 
2/N, the residue coming from the ground state energy 
is again the leading contribution to the partition func- 
tion as T ^ 0. If one calculates the magnetisation per 
particle, then, as conjectured, the decrease of M/N with 
increasing N at fixed T is suppressed. However, the ana- 
lytic result obtained is of limited use for large N, because 
the restriction on J means one is simultaneously taking 
J to zero. 

The second regime is the strong coupling limit J ~> oo. 
In this limit, the system is effectively quenched, so we 
expect to observe little impact from quantum entangle- 
ment. When we calculate the partition function, we find, 
for J/^B ^ 1, that the main contribution to Z is due to 
the residues arising from the 2n+l states that have s = n. 
Furthermore, the energy difference between states of dif- 
fering s can be treated as independent of /i_B. Specifically 
we have 
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It is now straightforward to determine the magnetic sus- 
ceptibility: 
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That this result is independent of N is not surprising 
since we are considering interacting spins on a complete 
graph in the strong coupling limit. The effect of quantum 
entanglement is thus washed out in this limit. 

In summary, we have obtained an expression for the 
partition function of a finite quantum system in thermal 
equilibrium, when the equilibrium states are obtained by 
maximising the Shannon entropy on the space of pure 
states. We applied this result to study the thermal expec- 
tation of the magnetisation and magnetic susceptibility 
of a system of TV spin-i particles. In the case of nonin- 
teracting spins, we were able to determine explicitly the 
effect arising from quantum entanglement. In reality, in- 
dividual particles forming magnetic substances are not 
likely to be fully entangled. Nevertheless, there are par- 
tial entanglements within a system, and entanglements 
appear to have the effect of enhancing phase transitions. 
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